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We study sudden quantum quenches in which the initial states are selected to be either eigenstates 
of an integrable Hamiltonian, which is non-mappable to a noninteracting one, or a nonintegrable 
Hamiltonian while the Hamiltonian after the quench is always integrable (and mappable to a nonin- 
teracting one). By studying weighted energy densities and entropies, we show that quenches starting 
from the nonintegrable (chaotic) eigenstates lead to an "ergodic" sampling of the eigenstates of the 
final Hamiltonian, while those starting from the integrable eigenstates do not (at least for the sys- 
tem sizes accessible to us). This goes in parallel with the fact that the distribution of conserved 
quantities in the initial states is thermal in the nonintegrable cases and nonthermal in the integrable 
ones, and means that, in general, thermalization occurs in integrable systems when the quench starts 
from an eigenstate of a nonintegrable Hamiltonian (away from the edges of the spectrum), while 
it fails (or requires much larger system sizes) for quenches starting at integrable points that are 
non-mappable to free models. We test those conclusions by studying the momentum distribution 
function of hard-core bosons after a quench. 
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I. INTRODUCTION 

The experimental realization of ultracold quantum 
gases in (quasi-)one-dimensional (ID) geometries has 
stimulated much theoretical research on ID systems [ljj. 
In equilibrium, they are particularly interesting because 
of the enhanced role played by quantum fluctuations. In 
addition, when taken out of equilibrium, they mav 
hibit unique behavior such as lack of thermalization 
We note that those systems are well isolated by the ul- 
tra high vacuum in which they are confined, i.e., their 
dynamics is to a very good approximation unitary 
Hence, by lack of thermalization we mean the fact that, 
after perturbing the system and letting it evolve, experi- 
mental observables relax to time-independent values that 
are not described by conventional statistical ensembles. 
This effect has been attributed to the proximity to inte- 
grable points H|-[8j], and has motivated an extensive ex- 
ploration of the dynamics of integrable systems [8H27j . 

The most common dynamics studied in this context are 
the ones that occur after the so-called sudden quenches. 
In a sudden quench, the initial state \^>i) is selected 
to be an eigenstate \<p n ) of an initial Hamiltonian Hi 
(Hi\<p n ) = e n \4> n )). Then, at time r = 0, some parame- 
ter (s) of Hi is (are) changed instantaneously Hi — > Hp, 
and the dynamics is followed under the final (time- 
independent) Hamiltonian Hp {Hp\ip a ) = E a \ip a )) 



can be written as 



|tf ( r )) = e - i&FT,h \^i) 
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where C a = (ipa^i) are the overlaps of the initial state 
with the eigenstates of the final Hamiltonian. 

Within this protocol, the time evolution of the expec- 
tation value of an observable 6, (6(r)) = (*(r)|0|*(r)), 



(O(r)) 



c a \ 2 o aa + £ c* a c^~ E ^^ h o, 



a/3, 



(2) 

where Q( g = (ij} a \0\%l)p). If the observable relaxes to 
a time-independent value (up to fluctuations that van- 
ish and revival times that diverge in the thermodynamic 
limit), that value can be computed (up to corrections 
that vanish in the thermodynamic limit) by taking the 
infinite-time average of Eq. (|2]l 



(O(t)) = Jim 



1 



dr(6(T)) 
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where we have assumed that there are no degeneracies, 
or that they can be neglected. The infinite-time average 
can be thought of as the prediction of an ensemble, the 
"diagonal ensemble" Q, where the density matrix is di- 
agonal in the basis of eigenstates of the final Hamiltonian, 
i.e., /5de = E„ \C a \ 2 \^ a )(ip a \- 

Several theoretical studies have shown that in inte- 
grable systems, and for observables of interest in current 
experiments and models, (O)de is different from the pre- 
dictions of conventional statistical mechanics ensembles 
UM El EE HI- Among the (few-body) observables 
that have been studied are the density and momentum 
distribution function of hard-core bosons , the dou- 

ble occupancy of fermions [l8j , and spin-spin correlations 
in spin chains [ID, EE] ■ If one wants to go beyond con- 
ventional statistical mechanics ensembles to take into ac- 
count that nontrivial conserved quantities {/•,} (whose 
number scales polynomially with system size) are present 
in integrable systems, and maximizes the entropy under 
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those constraints [28|, |29[ , then a generalized Gibbs en- 
semble (GGE) is obtained [|[. Remarkably, the GGE 
has been shown to provide the correct description for the 
previously mentioned observables (and others) after re- 
laxation following a sudden quench ||-E EE EE] ■ 
The (grand-canonical) density matrix for the GGE can 
be written as M 



1 



Pgge 



>: \ / 



(4) 



where ^gge = Tr[e~^j is the partition function, Ij 
are conserved quantities, and Xj are their corresponding 
Lagrange multipliers. The latter are computed using that 
(* / |J,|* 7 )=Tr[p GGE i J ]. 

The lack of thermalization of few-body observables of 
interest in integrable systems has been attributed to the 
failure of the eigenstate thermalization hypothesis 0- 
0, E, US HH] ■ That is to say that, within a microcanon- 
ical window, the matrix elements O a p in an integrable 
system fluctuate over a finite range of values. Hence, 
the outcome of the quench dynamics depends on which 
cigenstatcs of the final Hamiltonian have a finite over- 
lap with the initial state. If the initial state samples the 
cigenstates within the microcanonical window in an un- 
bias way, i.e., "ergodically" , then for that initial state the 
integrable system will thermalize. This does not happen 
for most quenches between inte grab le s yste ms studied in 
the literature @-Q3 EE E E, HHIl • HoweverMt 
has been shown to occur for special initial states [llj 
with particular entanglement properties [l6j . as well as 
for initial thermal states at infinite temperature fl2j . In 
the latter work, finite temperatures were shown to lead 
to a bias sampling of the final eigenstates of the Hamil- 
tonian after a quench. 

In this work, we study what happens in quenches from 
initial states that are eigenstates of an integrable Hamil- 
tonian that is non-mappable to a noninteracting one, 
while the final Hamiltonian is integrable (and mappable 
to a noninteracting one). This is of interest as most stud- 
ies in the literature have focused on quenches where both 
the initial and final Hamiltonians are mappable to non- 
interacting ones. We also study the case in which the 
initial state is an eigenstate of a nonintegrable Hamilto- 
nian and the final Hamiltonian is integrable. As argued 
in Ref. [32j . initial states that are eigenstates of nonin- 
tegrable Hamiltonians (away from the edges of the spec- 
trum) can lead to thermalization even though the final 
Hamiltonian is integrable. This can be understood as the 
two-body interactions that break intcgrability in models 
of current interest lead to quantum chaos even in the 
absence of randomness [33l436| . similar to what happens 
in systems with two-body random interactions [37H40I ) . 
As a result, the eigenstates of nonintegrable Hamiltoni- 
ans become random superpositions of eigenstates of the 
integrable ones to which the quench is performed, i.e., 
they provide the un-bias sampling required for thermal- 
ization after a quench to intcgrability. We do not find 



indications that this happens in quenches between the 
integrable systems studied here. 

The exposition is organized as follows. In Sec. [Hi 
we introduce the model Hamiltonian and describe the 
quenches to be studied. The observables and ensembles 
of interest are introduced in Sec. IIIII The remainder sec- 
tions are devoted to the discussion of the results for: (i) 
overlaps between initial states and eigenstates of the fi- 
nal Hamiltonians (Sec. IIV|) . (ii) entropies (Sec. [Vj), (iii) 
conserved quantities (Sec. IVij) . and (iv) momentum dis- 
tribution functions fSec. IVIlf . A summary of the results 
is presented in Sec. I Villi 



II. MODEL HAMILTONIAN AND QUENCHES 

We focus our study on ID lattice hard-core bosons in 
a box (a system with open boundary conditions) with 
nearest-neighbor (NN) hopping t and interaction V, and 
next-nearest-neighbor (NNN) hopping t' and interaction 
V. The Hamiltonian reads 



L-l 



6 = E + H - c -) +v(fk-l 



i=l 



L-2 



E-*'(^+2+H.c.) + ^ (fa -I 



1=1 



2 J V n * +2 -2'' (5) 



where b\(pi) is the hard-core boson creation (annihila- 
tion) operator, hi = b\bi the number operator, and L is 
the lattice size. In addition to the usual commutation 
relations = Sij, b\ and bi satisfy the constraints 

b\ 2 = tij = 0, which preclude multiple occupancy of the 
lattice sites. In what follows, the NN hopping parame- 
ter sets the energy scale (t = 1), we take h = &b = 1 
(where fcs is the Boltzman constant), and the number of 
particles N is always selected to be TV = L/3. 

For t' = V' = 0, Hamiltonian <j5j) is integrable (it can 
be mapped onto the spin- 1/2 XX Z chain) In ad- 
dition, if V = 0, the model can be exactly solved by 
first being mapped onto the spin-1/2 XY chain through 
the Holstein-Primakoff transformation [4l| and then onto 
noninteracting fcrmions by means of the Jordan- Wigner 
transformation [42| . Hence, the many-body eigenstates 
of the fermionic counterpart are Slater determinants. 
Properties of Slater determinants allow one to calculate 
off-diagonal correlations of hard-core bosons in polyno- 
mial time for both pure states p3l - l45| and the grand- 
canonical ensemble [4y]. The presence of nonzero NNN 
terms break the integrabilityof the model and drive the 
system to a chaotic regime |33l - t36| . 

Here, we consider two types of quenches. Quench type 

I is within integrable models, i.e, we keep t' = V = at 
all times, and select the initial state to be an eigenstate of 
Hamiltonian © with Vi > 0, while the evolution is stud- 
ied under Hamiltonian ([5]) with Vp — 0. Quench type 

II is from a nonintegrable model to an integrable one: 
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namely, the initial state is selected to be an eigenstate 
of Hamiltonian §5§ with t\ = V T ^0 while the dynamics 
is studied under Hamiltonian ([5|) with t' F = V' F = 0. In 
Quench type II, we keep V = at all times. Note that 
the final Hamiltonian in both types of quenches is the 
same. It can be fully diagonalized for large systems by 
being mapped onto nonintcracting fermions as mentioned 
before. In addition, since V and t' = V ^ are only 
present in the initial Hamiltonian, in what follows we 
drop the sub-index "7" from Vj and t\ = V{. 

In order to compare results of quenches from eigcn- 
states of different initial Hamiltonians, we select the ini- 
tial states \^i) such that after the quench the systems 
have an energy 

Ei = (tf/|if F |*/>, (6) 

which always corresponds to that of the canonical ensem- 
ble (CE) at a fixed temperature T G e (the temperature is 
taken to be Tqe ~ 2 throughout this work). This means 
that Ei = Ti{pceH f }, where p CE = e-" F / TcE /Z CE is 
the CE density matrix, and Z CE = Tr{e~ AF ^ TcE } is the 
partition function. The initial states, which are eigen- 
states of the t-V and t-t'-V' models, i.e., none of which 
is mappablc to a free model, are computed by means of 
the FILTLAN package that utilizes a polynomial filtered 
Lanczos procedure [471 ] . Since open boundary conditions 
are imposed, the only remaining symmetry at 1/3 filling 
and for V < 2, is parity. Here, we restrict our analyses to 
the even parity sectors, the dimension of which is about 
1/2 of the entire Hilbert space. The largest systems we 
study have L = 24, whose total Hilbert space dimension 
is D = ( 2 g 4 ) = 735 471. The even parity subspace for 
those systems has dimension D C vcn = 367 983. This is 
more than ten times larger than the largest Hamiltonian 
sector diagonalized in previous full exact diagonalization 
studies of similar systems @, 0, I331436} . In this work, the 
combination of utilizing Lanczos algorithm for comput- 
ing the initial state, and full exact diagonalization via the 
Bose-Fermi mapping for the final Hamiltonian, is what 
allows us to perform a complete analysis for larger Hilbert 
spaces. 

III. ENSEMBLES AND OBSERVABLES 

Given the fact that the initial states (l^/)) are pure 
states, and that the dynamics is unitary, the von Neu- 
mann entropy is zero at all times. In order to charac- 
terize our isolated systems, we make use of the diagonal 
entropy Sue- It is defined as the von Neumann entropy 
of the diagonal ensemble [48} 

5DE = -^|C Q | 2 ln(|C Q | 2 ), (7) 

OL 

and satisfies all thermodynamic properties expected of 
an entropy :48| . 5de has been recently shown to be con- 
sistent with the entropy of thermal ensembles for generic 
(nonintegrable) systems after relaxation (49} . 



Similar to the studies in Refs. [Ill, [l_2j, we compare the 
DE with the entropy predicted by the CE, the grand- 
canonical ensemble (GE), and the GGE. Namely, we also 
compute 

Ei 

S C e = lnZ CE + — — , (8) 

-tCE 

, Ei-fxN 

S GE = m Zqe H = , (9) 

S gge = lnZ GG E + J2 X j^i\Ij\^r), (10) 

3 

where Ei was defined in Eq. (j6|), and Z G e, Ice, N, and 
Z GGE were also defined before. In Eq. ©, the chemical 
potential and the temperature are computed such that 
Ei = Ti[p GE H F ] and N = Tv[p GE N], where p GE = 
e -(h f -iiN)/t ge / Zge jg thc GE density matrix, Z GE = 
T T { e -(HF-tiN)/T GE } is tne partition function, and N is 
the total number of particle operator. 

As mentioned before, all our quenches share the same 
final Hamiltonian (V = t' = V = 0), which is inte- 
grable. The conserved quantities {Ij} for the GGE in 
that model are taken to be the occupation number op- 
erators of the single particle eigenstates of the fermionic 
Hamiltonian to which hard-core bosons can be mapped. 
From this selection, it follows straightforwardly that the 
corresponding Lagrange multipliers can be computed as 
Xj = ln[(l - (tf/l^ltfi))/^/!/,-!*/}] g. Note that, by 
construction, the GGE has the same distribution of con- 
served quantities as the diagonal ensemble. We also cal- 
culate thc distribution of conserved quantities predicted 
by the GE, namely, (Ij)ge = Tr[/5 G E^]- 

In addition, we study the momentum distribution 
function of hard-core bosons in the DE ((rhk}r>E = 
Tr[p DE TOfc]); the GE ((m fe ) G E = Tr[p GE m fc ]), and thc 
GGE ((m fe )GGE = Tr[/5 GGE mfc]). fn k is thc diagonal 
part of the Fourier transform of the one-particle density 
matrix operator f>ij = o[bp 

1 L 

and it was previously studied in systems out of equilib- 
rium with the same Hp in Refs. (a-RC)}. 

IV. OVERLAPS 

As mentioned before, in integrable systems, the lack of 
eigenstate thcrmalization implies that the distribution of 
overlaps of the initial state with the eigenstates of final 
Hamiltonian \C a \ 2 plays a very important role in deter- 
mining the expectation value of observables after relax- 
ation [Eq. fl2J]. Therefore, we first study the values of 
\C a \ 2 as a function of E a in quenches type I and type 
II and compare them to the weights of the eigenstates of 
the final Hamiltonian in thc CE, e~ Ba / TcE /ZcE- 
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FIG. 1: (Color online) Coarse-grained plot of the weights 
in the DE [(a)-(c), (e)-(g)], and in the CE [(d) and (h)], 
|C Q | 2 and e _B ° //TcE /^cE, respectively, (a)-(c) Results for 
quench type I for different values of V . (e)-(g) Results for 
quench type II with different values of t' = V . (d) and 
(h) Weights in the CE, which correspond to effective tem- 
peratures Tce = 2.006 (Ej = -4.976) and Tce = 1.989 
(Ei — —5.011), obtained for the quenches with V = 1.28 and 
t' = V = 1.28, respectively. All results are for L = 24, N = 8, 
and the temperatures are Tce = 2.00 ± 0.01. The color scale 
indicates the number of states per unit area in the plot. 



In Fig. [1] we show density plots of the coarse-grained 
values of \C a \ 2 for different quenches type I [Figs. QJa)- 
[Tb)] and type II [Figs. me)-QIh)]. The corresponding 
results for the canonical weights (Tce ~ 2) are presented 
in Figs. QJd) and [Uh). One can clearly see in Fig. Q] 
that, for all values of V in quench type I and t' = V' in 
quench type II, the distribution of weights in the DE is 
qualitatively distinct from that in the CE, which shows 
a simple exponential decay. In the diagonal ensemble, 
the weights between eigenstates that are close in energy 
fluctuate wildly (note the logarithmic scale in the y axes) , 
and a maximum in the values of |C Q | 2 [particularly clear 
in panels (a), (b), (e) and (f)] can be seen around the 
energy of the initial state after the quench Ej w —5.0. 
Hence, in contrast to the canonical ensemble, eigenstates 
of Hf whose energies are close to Ej usually have the 
largest weights. This was also seen in Ref. |50 |. 

Another feature that is apparent in Fig. [T] is that while 
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FIG. 2: (Color online) Weighted energy density functions 
P(E) for both the DE and the CE. (a)-(c) Results for 
quenches type I. (e)-(f) Results for quenches type II. The 
dashed (blue) lines depict the energy shell, where (a) Ei = 
-4.99, SE = 0.159, (b) = -4.99, 5E = 0.571, (c) 
Ei = -4.98, 5E = 1.08, (d) Ei = -5.00, SE = 0.253, (e) 
Tj = -5.00, SE = 0.936, and (f) Tj = -5.01, SE = 2.66. 
Solid lines show the results of the fits of P(E) to the most 
appropriate functional form. Thin lines are fits to a Breit- 
Wigner function with (a) E = -5.01, T = 0.042, _(b) 
E = -5.04, T = 0.101, (c) E = -5.17, T = 0.567, (d) E = 
—4.99, F = 0.187, while thick lines are fits to a Gaussian with 
(e) E = -5.16, A = 0.650, and (f) E = -5.70, A = 2.31. In 
all cases L = 24 and N = 8. 



the weights |C a | 2 decrease rapidly as one moves away 
from Ej, the number of states with those weights in- 
creases dramatically for E a > Ej due to the increase of 
the density of states. For this reason, the appropriate 
quantity to quantify how different regions of the spec- 
trum contribute to the observables in the diagonal en- 
semble (and in the CE) is the weighted energy density 
function 



P(E) = ^ E W ai 

\E-E a \<AE/2 



(12) 



where W a is |C Q | 2 for the DE and e - E "/ TcE /Z C e for the 
CE. The summation is limited to the states in a small 
energy window \E — E a \ < AE/2, the width of which 
AE is selected to be the same for coarse-graining in Fig. 
[TJ We have checked that our results are robust for the 
selected values of AE. 

The weighted energy density functions for the quenches 
studied in Fig. [T] are shown in Fig. [21 There one 
can see that in all cases, except for the strongest type 



5 



II quenches, P(E) in the CE is much broader than 
that in the DE. For comparison, we have also plotted 
the results for the energy shells, which are Gaussians 
(V^SE)- 1 exp[-(E-Ei) 2 /(2SE 2 )] with the same mean 
energy Ej and energy width SE 2 = J2 a \C a \ 2 {E a — Ej) 2 
as the diagonal ensemble. The energy shell has the max- 
imal number of cigenstatcs of the final Hamiltonian that 
is accessible to the particular initial state selected. In 
quenches type I [Figs. [2ja)-[2jc)], even for the largest V, 
P{E) in the diagonal ensemble remains narrower than 
the energy shell. This is an indication of the lack of "er- 
godicity" of the initial states associated with quench type 
I. On the contrary, in quenches type II [Figs. [2jd)-[2jf)], 
the width of P{E) in the diagonal ensemble approaches 
that of the energy shell with increasing t' = V'. For the 
strongest quench considered, t' = V' = 1.28, P{E) in the 
DE and the CE become very close to each other as well as 
to the energy shell [Fig. [^f)]. The latter shows that ini- 
tial states that are eigenstates of nonintegrable Hamilto- 
nians sufficiently distant from an integrable point fill the 
energy shell ergodically after a quench to integrability. 

To better characterize P(E) as one goes from the weak- 
est to the strongest quenches, we have fitted it to two 
different functions: (i) a Brcit-Wigner function 



P{E) 



and (ii) a Gaussian 



1 



2tt (E - E) 2 + (r/2) 2 



P(E) = ■ exp 

V^rA 



-{E-Ef 



2A 2 



(13) 



(14) 



where the mean energy E and half-width V of the Brcit- 
Wigncr function, as well as the mean energy E and half- 
width A of the Gaussian, arc taken as fitting parameters. 

Figures [2ja)-[2jc) show that, as the strength of the in- 
teraction increases in quenches type I, the weighted en- 
ergy densities in the DE are better described by Brcit- 
Wigner functions. In quenches type II [Figs. [2Jd)-[2If)], 
on the other hand, weighted energy densities are de- 
scribed by Brcit-Wigner functions for weak quenches and 
transition to Gaussians as the strength of the quench in- 
creases. The fact that P{E) in the latter quenches ap- 
proach Gaussians (and ultimately the energy shell) hints 
the possibility of observing thermalization in those cases. 
This is consistent with the results in Refs. [IUdll, which 
considered a different set of quench protocols and were 
averaged over different initial states. Note that, in our 
results, no average has been introduced. 

Using the Breit-Wigner function to fit the weighted 
energy densities for the weakest quenches of type II, as 
well as the quenches of type I, is motivated by analytic re- 
sults obtained in the random two-body interaction model 
[37l - l40j . and by recent numerical results for models sim- 
ilar to those considered here [H, [36|. In Refs. [H, HHi 
the weighted energy densities (called strength functions) 
were computed for quenches that were the reverse from 
those considered here, namely, starting from a state |$) 
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FIG. 3: (Color online) P(e) for (a)-(c) type I quenches and 
(d)-(f) type II quenches. They were obtained by averaging 
over the nine even states \ipa) that are closest in energy to 
Ei. Before averaging, we shift es of each distribution to the 
one obtained for the state \ip a ) that is closest in energy to 
Ei. The dashed (blue) lines depict the energy shell, where 
(a) e s = -3.72, Se = 0.139, (b) e s = -3.66, Se = 0.557, (c) 
e s = -3.59, 5e = 1.11, (d) e s = -3.76, 5e = 0.221, (e) e s = 
-3.83, 5e = 0.885, and (f) e s = -4.10, Se = 3.54. Solid lines 
show the results of the fits of P(e) to the most appropriate 
functional form. Thin lines are fits to a Breit-Wigner function 
with (a) e = -3.70, T = 0.074, (b) e = -3.79, T = 0.126, (c) 
e = -3.72, r = 0.308, (d) e = -3.82, T = 0.94, while thick 
lines are fits to a Gaussian with (e) e = —3.95, A = 0.578, 
and (f) e = -4.37, A = 3.54. In all cases L = 18 and N = 6. 



that is an eigenstate \tp a ) of an integrable Hamiltonian 
(in one case, our Hp) and projecting that state onto the 
eigenstates |r/> ra ) of the t-V and t-t'-V' Hamiltonians (our 
Hi), i.e., |$) = J2 n c n\4>n)- In our notation, the weighted 
energy densities from those works can be written as 



P(e) 



1 



:„|<Ae/2 



(15) 



where the sum runs over states in the spectrum of Hi 
within a small window of energy Ae. Correspondingly, 
the energy shell associated with P(e) is now a Gaus- 



e n whose width is r5e 



sian centered at eg 
£Jc„| 2 (e„- £s ) 2 . 

For chaotic systems, it was shown in Refs. [35| - |4C)| . 
that as the interaction strength increases, and the aver- 
age strength of the coupling between unperturbed eigen- 
states becomes of the order of the average unperturbed 
level spacing, P(e) transitions from a Breit-Wigner func- 
tion to a Gaussian. This transition also seen numeri- 
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cally in integrable systems [3|1 [36| through an analysis 
of states whose energy was in the center of the spectrum, 
i.e., whose Tqe = oo. In order to make contact with 
the results in Refs. [HI, [Hj], we compute P(e) for our 
quenches and in the regions of the spectrum relevant to 
out work. Note that since full exact diagonalization is re- 
quired to obtain all eigcnstates \4> n ), the largest systems 
considered here in the calculations of P(e) have L = 18 
and N = 6. Finite size effects are strong for those lattice 
sizes, so we average the results of P(e) over the nine even 
parity states \ip a ) that are closest in energy to Ej 

In Fig. [3l we show results of P(e) for the same Hamil- 
tonian parameters for which we previously studied the 
weighted energy densities. Comparing Fig. [3] and Fig. 
[2j one can see that the results for P(e) and P(E) are 
qualitatively similar. For both types of quenches, P(e) 
becomes broader with increasing V (integrable) and t' = 
V' (nonintegrablc). However, when Hi is integrable, 
P(e) remains better described by Brcit-Wigncr functions, 
while a transition from Breit-Wigner to Gaussian behav- 
ior, as well as a filling of the energy shell [better seen in 
Fig- EJf ) than in Fig. |2{f)], is only observed for a nonin- 
tegrable i/j. 

Our results for integrable systems are in contrast to 
those reported in Refs. [3o1. l36l| . where a transition from 
Breit-Wigner to Gaussian was also observed for an in- 
teraction quench within integrable systems. This can be 
attributed to the combination of two effects. One is the 
fact that the states \ip a ) selected here are far from the 
middle of the spectrum, the case in Ref. [H, [3(| . Within 
the entire spectrum, the mean level spacing in the mid- 
dle is the minimal one, i.e., the same perturbation may 
couple more states in that region than away from it. If 
the average level spacing is the only reason for the differ- 
ence, then a transition should be seen in our case if one 
increases the system size. 

However, there may be another reason which would 
make the differences remain in the thermodynamic limit. 
The result of avera ging over eigenstates in the middle 
of the spectrum [35l l36l| leads to an ensemble at infinite 
temperature. In that ensemble, the distribution of the 
quantities that are conserved after the quench could be 
featureless. This would result in an un-bias sampling 
of the eigenstates of the Hamiltonian after the quench, 
which may not occur for quenches starting at a finite 
temperature even if they are strong quenches. We have 
found this to be the case for quenches between hard-core 
boson systems in which one starts from a thermal state 
Only the initial state at the infinite temperature 
was found to ergodically sample the eigenstates after the 
quench, and a finite size scaling analysis showed that 
this remains true with increasing system size. We will 
come back to this point in Sec. IVI1 where the conserved 
quantities are studied for all quenches considered here. 



V. ENTROPIES 

The fact that the weighted energy density appears to 
be Gaussian after a quench to an integrable model docs 
not automatically guarantee that thermalization will oc- 
cur, as it does not mean that an un-bias sampling has 
been performed (unless the energy shell is filled). One 
needs to keep in mind that a coarse-graining is involved 
when calculating P(E), where AE [see Eq. (fT2)) ] is much 
larger than the average level spacing. For example, in 
Ref. we showed that quenches between integrable 
systems, which started from thermal states, led to Gaus- 
sian like weighted energy distributions that are not ther- 
mal and do not become thermal with increasing system 
size. The exception was the initial infinite temperature 
state, which does lead to a thermal distribution after 
a quench (but at infinite temperature). A qualitatively 
similar behavior was seen for quenches starting from spe- 
cial pure states that are ground states of an integrable 
Hamiltonian (Tlj . 

In order to quantify how the weights evolve in the DE 
as the system size increases, and how they compare to 




V t'=V 

FIG. 4: (Color online) S/L for the DE, CE, GGE, and GE in 
systems with L — 15, 18, 21, and 24. Results are presented 
for (a)-(d) quenches type I and (e)-(h) quenches type II. 
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FIG. 5: (Color online) (Sge — Sde)/L as a function of (a) V 
for quench type I and (c) t' = V' for quench type II. (Sge — 
Sgge)/L as a function of (b) V for quench type I and (d) 
t' = V' for quench type II. Four different system sizes are 
presented for comparison. 



the weights in the GGE, CE, and GE, we have computed 
the entropies of the corresponding ensembles for all the 
quenches studied before. In Fig. |H we show all those 
entropies, per site, as a function of V for quenches type 
I [Figs. IHa)-[3Jd)] and of t' = V for quenches type II 
[Figs.0Te)-[l|h)], for four system sizes. Since, in all cases, 
we consider a fixed effective temperature Tqe ~ 2 after 
the quench, the values of S/L in the CE and the GE are 
(almost) independent of the quench parameters. 

The first feature that is apparent in Fig. @] is that, 
in both types of quenches, the entropy in the DE and 
the GGE are very small for weak quenches and increase 
with increasing quench strength. The GGE entropy fol- 
lows (but is always above, as expected from its grand- 
canonical character) the DE entropy. In addition, for 
each quench, the difference between Sde, Sgge and 
ScE, Sge is seen to decrease with increasing system size. 
However, only in quench type II one can see Sgge to 
become practically indistinguishable from Sge- If one 
agrees that the GGE describes observables after relax- 
ation, then the agreement between Sgge and Sge im- 
plies that the observables thermalize in those quenches. 

A better understanding of how the differences be- 
tween entropies scale with increasing system size can be 
gained through Fig. [SJ There we plot (Sge — Sue)/L 
and (Sge — Sqge)/L for different system sizes. Since 
(Sge — £>ce)/L vanishes in the thermodynamic limit, we 
drop Sge from the remaining analysis. We note that the 
analysis in Fig. [5] takes into account small fluctuations in 
Sge (unnoticeable in Figs. @| that occur because Ej is 
close but not identical in all initial states. 

Figure [5] shows that, with increasing system size, a 
decrease of (Sge — Se>e)/L and a vanishing of (Sge — 
Sqge)/L is apparent only in quenches type II, which, 
once again, indicates that thermalization will occur in 
those quenches. Furthermore, the value of t' = V at 



which an abrupt reduction of the differences (Sge — 
Sde)/L and (Sge — Sgge) /L occurs decreases as the sys- 
tem size increases, suggesting that in the thermodynamic 
limit an infinitesimally small type II quench will lead to 
thermalization. We should add that, in Fig. 03d), the 
slight upturn (note the logarithmic scale in the y axes) 
of (Sge — Sqge)/L for the strongest quenches seems to 
be related to the skewing of the weighted energy density 
seen in Fig. [2] Its effect is imperceptible in Fig. O^c), as 
(Sge — Sde) I L is more than an order of magnitude larger 
than (Sge — Sqge)/L. From the evolution of the data 
with increasing size, we expect this upturn to disappear 
in the thermodynamic limit. For quenches type I, the 
results in Figs. [Ha) andJSfb) show that (Sge — Sde)/£ 
and (Sge — Sgge)/L, respectively, either remain finite 
in the thermodynamic limit or vanish very slowly with 
increasing system size. 



VI. CONSERVED QUANTITIES 

As mentioned before, the presence of nontrivial sets 
of conserved quantities make integrable systems differ- 
ent from nonintegrable ones. Hence, an understanding of 
whether a particular initial state leads to thermalization 
can also be gained from analyzing how the quantities 
that are conserved after the quench behave in the ini- 
tial state. If the distribution of conserved quantities is 
identical to the one in thermal equilibrium with energy 
Ei after the quench, then the initial state can provide an 
un-bias sampling of the thermal one and lead to thermal- 
ization. Two particular examples in which that happens, 
in quenches between integrable systems, were discussed 
in Refs. [HII3- However, in those examples, the distri- 
butions of conserved quantities were featureless and cor- 
responded to systems that were at infinite temperature 
after the quench. 

In Fig. |5l we show the distribution of conserved 
quantities for various initial states for quenches type I 
[Figs. (Bfa)-^d)] and type II [Figs. [HIe)-|n][h)], and com- 
pare them to the distribution of conserved quantities in 
the GE (which is almost identical in all quenches as Ei 
is very close in all of them). That figure shows that 
in quenches type I remains different from that 
in the CE for all values of V. In contrast, in type II 
quenches, (ty i) becomes very similar to that in the 
CE as t' = V increase. 

In order to have a more quantitative understanding of 
the difference between the conserved quantities in the ini- 
tial state (in the GGE) and in the GE, as the system size 
increases, we calculate the relative integrated difference 
defined as 



AI = 



Kgf&ljEj) - (^)ge| 



(16) 



Results for AI in both types of quenches, and for dif- 
ferent systems sizes, are presented in Fig. [5] They make 
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FIG. 6: (Color online) Distribution of conserved quantities 
in the initial state (same as the GGE) and in the 
GE. Results are presented for (a)-(d) type I quenches, and 
(e)-(h) type II quenches. In all cases, L = 24 and N — 8. 



evident that, as the system size increases in quenches 
type II, the conserved quantities in the initial state con- 
verge to those predicted in thermal equilibrium. No such 
clear tendency is seen for quenches type I. We note that 
the results for AI are qualitatively similar to those ob- 
tained for (Sge — Sqge)/L in Fig. [7] This can be un- 
derstood as the entropy in the initial state (GGE) and in 
the GE can be written in terms of the occupation of the 
single-particle eigenstates in each case 



S 



L 

■Eta 



/,i /J"//. (17) 



where Ij = for the initial state (GGE) and 

Ij = (4>GB for the GE. 

The results obtained for quenches type II are remark- 
able as they imply that such quenches lead to thermal- 
ization in integrable systems. This occurs even though 
the distribution of conserved quantities is a nontrivial 
one (it is not flat, which was the case in Refs. (111. fl2j). 
Beyond its implication for the quantum dynamics of inte- 
grable systems, one could think of using this information 
to learn about complicated many-body systems. It im- 
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FIG. 7: (Color online) Relative integrated difference of the 
conserved quantities in the initial state (in the GGE) and the 
GE, AI, as a function of (a) V for quench type I and (b) 
if = V' for quench type II. Results are presented for four 
system sizes. 



plies, e.g., that if we know the kinetic energy of a chaotic 
strongly correlated fermionic system in equilibrium, we 
could automatically calculate its momentum distribution 
by computing the momentum distribution function of a 
system that is in thermal equilibrium in the noninteract- 
ing limit with the same energy as the kinetic energy of 
the strongly correlated one. 



VII. MOMENTUM DISTRIBUTIONS 

In this section, we test whether the conclusion reached 
in the previous sections, that quenches type I do exhibit a 
clear tendency to thermalize with increasing system size 
while quenches type II do, holds for an observable, the 
momentum distribution function rrifc. 

In Fig. [51 we show the momentum distribution func- 
tion for the diagonal ensemble in various initial states for 
quenches type I [Figs. |5Ja)-|5I[d)] and type II [Figs. |5Je)- 
|5Jh)], and compare them to the momentum distribution 
functions predicted by the GGE and the GE. That fig- 
ure shows that while the GGE results closely follow the 
ones in the DE for all quenches, the CE results are only 
consistently closer to the DE ones in quenches type II as 
one increases t' = V [Figs. |5{f)-(Sfh)] • 

Once again, in order to be more quantitative, we com- 
pute the integrated differences between the GE and DE 
results, namely 



AfflDE 



(rhk) 



DE 



Efe(w/c)GE 



(18) 



Results for Atode are presented in Fig. Eta) for 
quenches type I and in Fig. EJb) for quenches type II. In 
the former figure, no consistent trend is seen in the Atode 
with increasing system size. In contrast, in Fig. |9]Jb) one 
can see that, for the three largest system sizes, Atode 
decreases steadily with increasing L for t' = V between 
0.1 and 1. This is in agreement with what was expected 
from our previous analyses, and indicates that quenches 
type II lead to thermalization for this observable, while 
quenches type I do not lead to thermalization or require 
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FIG. 8: (Color online) Momentum distribution functions in 
the DE, GGE, and GE. Results are presented for (a)-(d) type 
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FIG. 9: (Color online) Relative integrated difference Attide 
as a function of (a) V for quench type I and (b) t' — V' for 
quench type II. Results are presented for four system sizes. 



state is selected to be either an eigenstate of an inte- 
grable Hamiltonian that is non-mappable to a free model 
(quench type I) or an eigenstate of a nonintegrable one 
(quench type II). By studying weighted energy densi- 
ties and entropies, we have not found clear evidence that 
quenches type I, at least within our Hamiltonians of in- 
terest and away from the middle of the spectrum, lead to 
an un-bias sampling of the eigenstates of the final Hamil- 
tonian. Quenches type II, on the other hand, are found 
to provide such an un-bias sampling. Furthermore, an 
analysis of different systems sizes indicates that, in the 
thermodynamic limit, an infinitesimal quench type II will 
lead to thermalization. Here, an important requirement 
to keep in mind is that the initial state must be away 
from the edges of the spectrum. This is because, for sys- 
tems in the absence of randomness and with two-body 
interaction, chaotic eigenstates can only be found away 
from the edges of the spectrum [H, H|| . 

We have also shown that, in the initial state, an analy- 
sis of the distribution of the quantities that are conserved 
after the quench provides an understanding of why ther- 
malization occurs in one type of quenches while it fails 
in the other one. In quenches type I, that distribution 
remains different from, or approaches very slowly with 
increasing systems size, the one in thermal equilibrium 
after the quench. This implies that an un-bias sampling 
of the eigenstates of the final Hamiltonian does not occur 
or takes very large systems to be discerned. The oppo- 
site is true for quenches type II. Hence, quenches type 
II provide a consistent way of creating initial states that 
have the right distribution of conserved quantities so that 
thermalization can occur after a quench to intcgrability. 
Special initial states for which this occurred in quenches 
type I were discussed in Refs. [ill, EH- However, the dis- 
tribution of conserved quantities in those cases was (triv- 
ially) fiat corresponding to infinite temperature systems 
after the quench. 

Finally, by studying the momentum distribution func- 
tion in quenches type I and type II, we have shown that 
the conclusions reached on the basis of the results for 
the energy distributions, entropies, and conserved quan- 
tities hold. Namely, we have found no indications that 
quenches type I lead to the thermalization of this observ- 
able while quenches type II do result in thermal behavior. 



much larger system sizes to observe an approach to the 
thermal predictions. 



VIII. SUMMARY 

We have studied quenches to an integrable Hamilto- 
nian (mappable to free fermions) in which the initial 
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